Comparison of the Total Charged-Particle Multiplicity in High-Energy Heavy Ion 

Collisions with e^e~ and pp/pp Data 
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The PHOBOS experiment at RHIC has measured the total multiplicity of primary charged par- 
ticles as a function of collision centrality in Au-|-Au collisions at y^sjvjv ~ 19.6, 130 and 200 GeV. 
Above ^/s ~ 20 GeV, the total multiplicity per participating nucleon pair {{Nch) /{Npart/2)) in cen- 
tral events scales with in the same way as (Nch) in e^e~ data. This is suggestive of a universal 
mechanism of particle production in strongly-interacting systems, controlled mainly by the amount 
of energy available for particle production (per participant pair for heavy ion collisions). The same 
effect has been observed in pp/pp data after correcting for the energy taken away by leading par- 
ticles. An approximate independence of (Nch) /{Npart/2) on the number of participating nucleons 
is also observed, reminiscent of "wounded nucleon" scaling {Nch <x Npart), but with the constant of 
proportionality set by the multiplicity measured in e^e~ data rather than by pp/pp data. 

PACS numbers: 25.75.Dw 
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Central collisions of two gold nuclei at the top en- 
ergy of the Relativistic Heavy Ion Collider (RHIC) at 
Brookhaven National Laboratory produce thousands of 
charged particles. These are the largest particle mul- 
tiplicities generated in man-made subatomic reactions. 
The hope is that these complex systems may reveal evi- 
dence of the creation and decay of a Quark-Gluon Plasma 
(QGP), where quarks and gluons are allowed to explore 
a volume larger than that of a typical hadron. 

The high multiplicities in heavy ion collisions typically 
arise from the large number of nucleon-nucleon collisions 
which occur, with many of the nucleons struck several 
times as they pass through the oncoming nucleus. Stud- 
ies of proton-nucleus collisions demonstrated that the to- 
tal multiplicity (Nch ) does not scale proportionally to the 
number of binary collisions (Ncoii) in the reaction, but 
rather was found to scale more closely with the num- 
ber of "wounded nucleons" which participate inelasti- 
cally (Npart) For example, the number of par- 

ticipants is Npart = 2 for a proton-proton collision and 
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1) for a proton-nucleus collision. Thus, 



by scaling the particle yields measured in heavy ion col- 



lisions by Npart /'^^ data from heavy ion collisions may be 
directly compared with similar yields in elementary pp, 
pp or even the annihilation of e+e^ into hadrons. 

While both e'^e~ and pp/pp collisions must ultimately 
allow a description based on Quantum Chromodynamics 
(QCD), the theory of the strong interaction, the evo- 
lution of these two systems tends to be understood in 
different ways. The large momentum transfer to the out- 
going produced quark and anti-quark in e+e^ reactions 
allows the use of perturbative QCD (pQCD) to describe 
the spectrum of quarks and gluons radiated as the system 
fragments Q . Minimum bias collisions of hadrons are not 
generally thought to be amenable to such a perturbative 
description, since the transverse momentum exchanges 
involved are typically less than 1 GeV/c. Instead, phe- 
nomenological approaches (e.g. PYTHIA ^) are used to 
describe most of the (predominantly soft) particles pro- 
duced in high energy pp or pp collisions. 

A basic connection between perturbative and non- 
perturbative physics has been elucidated by simultaneous 
measurements of the multiplicity of charged particles and 
the high- momentum "leading" protons in pp collisions at 
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FIG. 1: (a) dNch/drj/ (Npartfi) of charged particles produced 
in central Au+Au collisions at ^/snn = 200 and 19.6 GeV 
(shifted by Ari — 2.32), compared with elementary systems. 
A Woods-Saxon fit to the 200 GeV Au+Au data is shown. 
The e'^e~ data are plotted as a function of yr, the rapid- 
ity relative to the thrust axis, (b) PHOBOS and UA5 data 
divided by a Woods-Saxon fit to the 200 GeV Au-I-Au data. 



the ISR. Basile et al. 5] found that the average multi- 
phcity (Nch) in PP collisions is similar to that for e+e^ 
collisions with y/s^+^- = ^/ScS, where ^yseS is the pp 
center-of-mass energy minus the energy of the leading 
particles. This is interpreted as a universal mechanism 
of particle production controlled dominantly by the avail- 
able center of mass energy 

In this Letter, we report results from the PHOBOS 
experiment on the total multiplicity of primary charged 
particles (Nch) as a function of Npart for heavy ion col- 



lisions at y/sNN — 19.6, 130 and 200 GeV, where ^/snn 
is the nucleon-nucleon center-of-mass energy. Compar- 
isons with pp/pp and e~^e~ data are made to investigate 
whether this universal mechanism of particle production 
applies in the context of heavy ion collisions. 

The PHOBOS multiplicity detector consists of two ar- 
rays of silicon detectors which cover nearly the full solid 
angle for collision events. The "Octagon" detector sur- 
rounds the interaction region with a roughly cylindrical 
geometry covering \ri\ < 3.2. Two sets of three "Ring" 
detectors are placed far forward and backward of the in- 
teraction point and surround the beam pipe, covering 
3 < jryl < 5.4. The methods used for measuring the 
multiplicity of charged particles as well as for extracting 
(Npart) have been described in more detail in Refs. 



Using the data presented in Ref. Fig. QJi shows 
dNch/dri/ (Npart) /2 averaged over the forward and back- 
ward hemispheres for the most central Au-|-Au events at 
y^SNN = 200 GeV. The systematic errors (representing a 
90% CL interval) depend on rj and are shown on the figure 
as a shaded band. The Au-fAu data are compared with 
dNch/drj for non-single diffractive (NSD) pp collisions Q 
and dN/dyT for e'^e~ collisions (with cuts applied to re- 
ject large initial-state photon radiation) at y/s = 200 
GeV. The variable yr is the rapidity of charged particles 
relative to the event thrust axis, assuming all particles to 
have the pion mass. 

It is observed that the Au+Au data are very similar in 
magnitude and shape to the e^e~ data at the same ^/s, 
and similar in shape to the pp data (as shown in Fig. ^p), 
over a large range in ij. The differences between the e'^e~ 
and Au+Au distributions shown in Fig. ^ can be partly 
attributed to the different kinematic variables. JETSET 
calculations indicate that the yx distribution is slightly 
narrower than the corresponding pseudorapidity distri- 
bution in e"'"e~ collisions, with a higher plateau height. 
Yet, even without taking this into account, the difference 
between the distributions is no more than ±10% for \r]\ 
and < 4 4j. However, the same calculations also 
show that the choice in kinematic variables does not ex- 
plain the difference in the forward region (above \ri\ — 4), 
although this may not be surprising, as this region should 
show some residual effect of the presence of the spectator 
nucleons. 

The similarity of the angular distributions indicates 
that the total yield of charged particles in e+e^ and 
central Au+Au collisions should also be similar for the 
same ^/s, when the nuclear data are scaled by the num- 
ber of participant pairs. To correct for the small ac- 
ceptance losses in the PHOBOS apparatus (which covers 
\ri\ < 5.4), we have used several methods inspired by the 
excellent agreement of the lowest energy PHOBOS data 
with the higher energy data when shown as a function of 
1]' = rj — ybeam PHOBOS data from ^snn — 19.6 
GeV for rj > 2.5, shifted by Ar] — y2m ^ yws — 

2.32 

(the difference in beam rapidities between the two ener- 
gies) , displays the limiting behavior discussed in Ref. . 
This effectively extends the rapidity coverage to ~ 8. 
A Woods-Saxon function for dN/dy fit to the Au+Au 
data, also provides a reasonable description of the dN/ drj 
distribution, and extrapolates through the lower energy 
data as well. Thus, in one method, we integrate dNch/drj 
for y/sNN = 130 and 200 GeV for 77' < and then use 
the PHOBOS data at Vsjvw = 19.6 GeV for 77' > 0. We 
also integrate Woods-Saxon fits, similar to that shown in 
Fig. for 1 77 1 < 8. These two approaches agree within 
2% for central events. For the lowest RHIC energy, we 
simply integrate the charged particles in the PHOBOS 
acceptance. 

In Fig. 13^, data on Nch from pp, pp, e+e~ and central 
heavy ion collisions (scaled by Npart f'^) are shown as a 
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function of y/s. The pp, pp, and e+e^ data and errors are 
taken from a compilation 10] and no further corrections 
are apphed. The errors shown arc the quadraticaUy com- 
bined statistical and systematic errors. Heavy ion data 
are shown for central Au+ Au events at RHIC (this work) , 
Au+Au events from E895 at the AGS {./s^n = 2.6-4.3 
GeV) [ni and Pb+Pb events from NA49 at the SPS 
(^/i^ = 8.6, 12.2 and 17.3 GeV) A PHOBOS 

Au+Au data point at ^smn = 56 GeV has been added 
by using the measured value at midrapidity [l^ and us- 
ing the universal limiting distribution described in Ref. 
[3| to approximate the shape of the full distribution. All 
of the errors shown for the heavy ion data are systematic. 

Perturbativc QCD calculations are able to predict the 
dependence of the total multiplicity in e^e~ collisions as 
a function of ^/s, iVe+e-(s) — (7q!s(s)"^ exp(-yB7a7(s)), 
with A and B fully calculable within pQCDui]. The 
QCD scale Aqcd is set to 225 MeV, leaving only a con- 
stant of proportionality C free to fit to the experimental 
data. A fit to the e+e" data has been made with this 
expression ("e+e~ fit") and has been used in Fig. |2h to 
see how the various systems compare with e^e~ data by 
scaling all of the data at a given ^/s by this function. 

Fig. 1213 shows that the pp/pp data are about 30% below 
e~^e~ over the full range of energies. However, rescaling 
the y/s of each point by a factor of 1/2, ^Soff = •\/s/2, 
brings the data into reasonable agreement with the e'^e^ 
trend, as shown by the open diamonds. This is consistent 
with measurements of leading protons in pp collisions, 
which find dN/dxp (where xp = '2.pz/\fs in the collider 
reference frame) to be approximately constant for non- 
diffractive events over a large range of -^/s [l5j | and thus 
{xf) - 1/2. 

Unlike the pp/pp data, the heavy ion data does not 
follow the e^e~ trend over the whole energy range. In- 
stead, they lie below the pp data at AGS energies, crosses 
through the pp data between AGS and SPS energies, and 
joins smoothly with the e^e~ data above the top SPS en- 
ergy. Thus, at high energies, the multiplicity measured 
per participant pair in Au+Au collisions scales in a sim- 
ilar way to e^e~ data at the same ^/s. If we understand 
the lower effective ^/s in pp collisions as stemming from 
the "leading particle effect", where the leading proton 
carries off a substantial amount of the available energy, 
the Au+Au data suggest a substantially reduced leading 
particle effect in central collisions of heavy nuclei at high 
energy. 

The alleviation of the leading particle effect might not 
be so surprising in central nuclear collisions. In the 
Glauber model, each participating nucleon is typically 
struck 4 — 6 times on average as it passes through the on- 
coming gold nucleus in a central event (the exact value 
of V depending on the energy-dependent nucleon- nucleon 
inelastic cross section, (Jnn{s)). One could speculate 
that the multiple collisions simultaneously excite and 
dissociate the participating nucleons, transferring much 
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FIG. 2: (a) The total charged multiplicity (Nch) for pp, pp, 
e^e~ and central Au+Au events as a function of y^. The 
Au+Au data are normalized by Npart/2. The dotted line is 
a perturbative QCD expression fit to the e'*'e~ data. The 
diamonds are the pp/pp data with ^Scft = \/s/2- (b) The 
data in (a) divided by the e^e~ fit, to allow direct comparison 
of different data at the same ^/s. 



more of the energy from the forward direction towards 
midrapidity than found in an average pp/pp collision - but 
ultimately limited by the total incident energy. This hy- 
pothesis should be testable in proton-nucleus collisions, 
by measuring particle yields as a function of v as was 
done in Refs. [lallS- The data in those references sug- 
gest that pion yields, whether in the projectile region 
{y > 0) integrated over 47r 0| increase rapidly 

for v < 3 and then much more slowly for v > 3. How- 
ever, limited experimental acceptances and theoretical 
uncertainties preclude making any strong conclusions re- 
garding the relationship between the energy loss of the 
projectile and the total charged multiplicity. 

In Fig.^iNch)/ {Npart/2) is shown for PHOBOS data 



at three RHIC energies as a function of Nj 



part • 



The 



90% CL systematic error on the centrality dependence 
of (Nch) / {Npart/2) is shown as a shaded band, and rep- 
resents a combination of several factors, dominated by 
the uncertainty of the extrapolation procedure to extract 
Nch over the full solid angle. 

It might have been expected that, in events with larger 
impact parameters, each participant would have fewer 
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collisions on average and thus not be fully dissociated. 
However, within the systematic errors, the total yield 
per participant pair is approximately constant (within 
10%) over the measured centrality range, 65 < (Npart) < 
358, which corresponds to 3 < F < 6, where V is the 
average number of collisions undergone by each oncoming 
nucleon. Thus, it appears that only the first few collisions 
have any appreciable effect on particle production. It 
should be noted that this simple scaling is not observed 
for particle yields measured in a limited pseudorapidity 
range near midrapidity 

Proton-antiproton data exist at 200 GeV, but not for 
the other two RHIC energies. We use a parametrization 
ofpp data from Ref. Hi, (N^h) = -4.2 + AMs '^^'^, for 
19.6 and 130 GeV. Several measurements exist in e"*"e~ 
at 200 GeV, but not for the other two energies. For 
these we use the pQCD formula for iVg+g- , the quality 
of the fit clearly indicated in Fig. ^p. Fig. |21 shows that 
the Au+Au data are consistent with "wounded nucleon" 
scaling, in that the multiplicity is proportional to Npart 
(Nch oc Npart)- However, Nch clearly does not scale sim- 
ply with the multiplicity measured in pp collisions at the 
same energy. Rather, for a large range of impact pa- 
rameter, the multiplicity scales approximately with the 
total multiplicity in e+e~ annihilation at the same ^/s. 
Thus, it appears that the first few collisions per partici- 
pant are sufficient to liberate as much energy for particle 
production as an e+e^ reaction. 

However, the rapid approach of (Nch) / {Npart/ 2) in 
central heavy ion collisions below ^/JWn ^ 20 GeV to- 
wards the e^e^ data complicates any simple geometric 
interpretation, as all of the heavy ion data compared are 
for a similar range of impact parameters. One feature 
that might point to why the particle yields at the AGS 
and SPS are perhaps "suppressed" relative to e'^e~ data 
(and even to pp data at lower energies, as noted in Ref. 
[131) is the magnitude of the ratio of net baryons to pions 
in the system. This ratio, which scales approximately 
as Npart/Nch, is 0(50%) at AGS energies 11], but is 
(1%) at RHIC [1^. In a thermal statistical approach 
this refiects the decrease of the baryon chemical po- 
tential with increasing beam energy. 

In conclusion, the PHOBOS experiment has mea- 
sured the normalized charged-particle multiplicity 
{Nch) I {Npart I "2) in Au+Au collisions at three RHIC en- 
ergies as a function of the centrality of the collision. 
Above CERN SPS energies, the total multiplicity per 
participating nucleon pair, {Nch) / {Npart/ '2) i in central 
events scales with -^/s in the same way as e+e~ data. 
This is suggestive of a universal mechanism of parti- 
cle production in strongly-interacting systems, controlled 
mainly by the amount of energy available for particle 
production. This may be related to the multiple colli- 
sions suffered by each participant nucleon, which could 
substantially reduce the leading particle effect seen pre- 
viously in pp collisions. The weak centrality dependence 
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FIG. 3: {Nch) /{Npart/ 2) is shown vs. Npart for ^sjvjv = 
19.6, 130, and 200 GeV as closed symbols. The error includes 
contributions from the uncertainty on overall Nch scale and 
Npart scale. The shaded band shows the uncertainty on the 
extrapolation procedure. The open symbols show UA5 data 
at 200 GeV and results from an interpolation for the lower 
energies. The dotted lines show the values from the e^e~ fit. 



for {Nch) /{Npart/ 2), reminiscent of "wounded nucleon" 
scaling, suggests that after the first few collisions per par- 
ticipant, the multiplicity per participant pair saturates 
near the value measured in e+e~ reactions. Ultimately, 
the existence of simple scaling behavior with ^Soff and 
Npart indicates stronger constraints on particle produc- 
tion than previously considered theoretically. Thus, these 
results may provide a new perspective on particle produc- 
tion in heavy ion collisions. 
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